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ABSTRACT 


A computer simulation of an aircraft making a low level penetra- 
tion over hostile terrain is described. The model discussed is intended 
as a tool for the determination of appropriate defensive tactics for air- 
craft upon receiving fire from conventional ground weapons. In any com- 
puter model of a ground to air combat environment it is essential to have 
a realistic simulation of the terrain. The method used in this simulation 


is discussed and shown to be extremely accurate. 
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i Introduction 


PENAIR is a time step computer simulation of an aircraft flight over 
hostile terrain. It is an analytic computer war game in which one aircraft 
and a maximum of two hundred weapons are the principles. Any piece of 
actual terrain may be simulated using a contour map of the area to obtain 


the input parameters. 


The simulation is intended as a tool for the determination of appropri- 
ate defensive tactics for an aircraft aan receiving fire from conventional 
ground weapons. Using the characteristics of the aircraft and of the wea- 
pons the model determines the survival probability of the aircraft on any 
of several flightpaths. There are no restrictions on the aircraft's speed, 
altitude or other characteristics and therefore this model is equally suita- 
ble for helicopters or jets. 

The aircraft flightpath can be predetermined or it can be generated 
by the model. Upon receiving fire the aircraft can be programmed to exe- 
cute any of 18 evasive maneuvers. A comparison of the survival probabili- 
ties can then be made thus giving an indication of the appropriate defensive 
maneuver to employ. 

The terrain simulation and generated flightpaths appear to be ex- 
tremely realistic. The terrain is simulated by fitting, in the least squares 
sense, a polynomial function of two independent variables (x,y) to values 
of a dependent variable specified at points on a rectangular grid in the 
plane of the independent variables. This is accomplished by a separate 
computer program, TERRAIN, which generates the necessary inputs for 


PENAIR. 


ZL. Design of the Model 
PENAIR is a computer simulation of a single aircraft making a low 
level penetration over terrain containing hostile small arms and con- 


ventional anti-aircraft weapons. Any actual or fictitious terrain may be 





simulated. No offensive aircraft weapons are considered although 18 dif- 
ferent defensive maneuvers Such as increasing speed and diving for the 
deck may be performed. Ten types of ground weapons may be used with 

a maximum of 20 weapons per type. The program computes the aircraft's 
time in view of each weapon, the probability of the aircraft surviving each 
weapon and the aircraft's total survival probability. The simulation runs 

in discrete time steps using any time interval desired. The program does 
not utilize Monte Carlo methods; aircraft survival probability for a speci- 
fied flightpath over a specified terrain is computed in a single run. Flow 
chart 1.0 of Appendix I, page 30 , indicates the general logic of the PENAIR 


program. 


PENAIR is written in FORTRAN-63 for the CDC 1604 computer. The 
PENAIR program is included as Appendix II. A separate program, TERRAIN, 
is utilized to supply inputs to PENAIR for terrain simulation. The terrain 
input program supplies the necessary data on punched cards for the PENAIR 
program. The TERRAIN program is also programmed in FORTRAN-63 for the 
CDC 1604 and is included as Appendix IV. 

The terrain used in the PENAIR program may be divided into as many 
as 12 separate squares for simulation. Each of the 12 individual squares 
is approximated by a separate polynomial. The 12 squares may be of any 
common size. Figure 1 illustrates the arrangement of these squares. The 
coordinate system has been established such that the lower left corner is 
the origin with increasing x to the right and increasing y upward. It is 
possible to use actual military grid coordinates since the origin may 
take any desired x and y values. 

A special PENAIR output is a set of punched cards to be used as 
input to a graph plot routine for an off-line plotter. This routine then draws 
the actual flightpath over the simulated terrain. Output from this routine 


is included in Section 6. 








FIGURE 1 


THE ARRANGEMENT OF THE INDIVIDUAL 
SQUARES OF TERRAIN SIMULATED IN 
THE PENAIR PROGRAM 


3. Aircraft Movement 

The PENAIR model has provisions for three types of simulated aircraft 
flightpaths. The first of these flightpaths allow the PENAIR program to 
select one of 18 different evasive maneuvers for the aircraft upon receiv- 
ing fire. The 18 evasive maneuvers consist of six basic maneuvers with 
three possible speed changes for each. The six basic maneuvers are: 
(1) continue flying straight ahead at the same altitude, (2) commence a 
climb straight ahead at the aircraft's maximum rate of ascent, (3) tum 
away from the weapon remaining at the same altitude, (4) turn away and 
commence a maximum rate climb, (5) dive for the deck and commence nap 
of the earth (contour) flight straight ahead, or (6) turn away from the 
weapon, dive for the deck and commence nap of the earth flight. For each 
of the above maneuvers the aircraft may decrease speed, remain at the 
same speed or increase speed. All of the above maneuvers are gener- 
ated using the input characteristics of the aircraft being simulated. It 


should be realized by the PENAIR user however that only one basic maneuver 


may be employed in one computer run. The maneuver employed is deter- 
mined by the numerical value of the input flag IMARK. 

The other two simulated aircraft flightpaths are preplanned. In the 
first of these, the flightpath is completely determined in advance and 
read into the PENAIR program from data cards. In the second, the x and y 
coordinates of points along the flightpath are read in to the PENAIR pro- 
gram. From these points the PENAIR program computes a nap of the earth 
flightpath for the aircraft over the simulated terrain. In both of these 
cases, the PENAIR program makes no deviation from the pre-planned flight- 
path and the aircraft makes no use of the evasive maneuvers. A complete 
description of the inputs needed to regulate aircraft movement is included 
in Appendix III. 

4. Determination of the Aircraft's Survival Probability. 

The major output of the PENAIR program is the probability of aircraft 
survival. This survival probability is computed as a function of weapnn 
type, range and time in view. The necessary inputs are weapon character- 
istics and the single time-interval kill probabilities for each weapon. 
Given these kill probabilities, the probability of an aircraft surviving any 


single weapon is computed as: 


P (Survival) = (1 - P,) (1- P.) anal a) 


where P. is the probability of kill for one time interval at rangei fora 
given weapon. The total probability of aircraft survival is then computed 
as the product of the individual probabilities of survival for each weapon. 
(See Flow Chart 5.1, Appendix I, page 48). It is assumed that the proba- 
bility of a kill by one weapon is independent of all other weapons. The 
number of time intervals the aircraft receives fire depends on the time in 
view, minimum and maximum weapon range, and the weapon acquisition, 


firing, and reloading times. 





Time in view is the one parameter to be determined. The method used 
by this simulation to determine time in view will be explained using Figure 
2 and Flow Chart 7.0 of Appendix I, page 49, 

The aircraft is considered in view of a weapon if every point on the 
line of sight WA between the weapon position W, and the aircraft A, is 
higher than the corresponding point of the terrain. To determine this the 
x,y coordinates of the point on the line TR, 10 meters from the weapon are 
used to calculate the elevation of the terrain at that point. This elevation 
is then compared with the height of the line of sight at those coordinates. 
If the terrain is lower, another point 50 meters farther along the line TR 
is chosen and the same comparison is made. If the terrain elevation is 
below the line of sight elevation at this point this procedure is repeated 
every 50 meters. This process is continued until it is determined that the 
aircraft is in view of the weapon or that the terrain elevation exceeds the 
line of sight. If the latter occurs the aircraft is considered not in view 
of this weapon at this time. This process is repeated for all weapons and 
then the aircraft is advanced one time interval. The number of consecutive 
time intervals the aircraft remains in view of any one weapon is then used 


to compute the current total amount of time in view per weapon. 





FIGURE 2 


DETERMINATION _ OF THE AIRCRAFT'S TIME 
IN VIEW OF A WEAPON 





Two other methods of determining if the aircraft was in view of the 
weapon were tested. Elevation comparisons between line of sight and 
terrain were made at 25 meter increments on the line of sight WA, and at 
25 meter increments along TR. No significant differences were noted in 
the times of view computed by the 3 methods mentioned. However, the 
method outlined in the previous paragraph greatly reduces computer rmn- 
ning time for the PENAIR program. 

After the aircraft has been in view of a weapon for the required acqui- 
sition time, the weapon minimum and maximum ranges are checked. If 
the aircraft is within firing range, the weapon commences firing and the 
probability of survival is determined for this weapon in this time interval. 
The probabilities of aircraft survival for each weapon are then computed 
for all following time intervals until the aircraft either goes out of view 
and/or range or until the weapon has expended its ammunition. In both 
cases no further probabilities are computed until the weapon is reloaded 
or until the aircraft comes back in view and within range. It is assumed 
that after reloading reacquisition requires one-half the acquisition time. 
If the aircraft goes out of view the weapon must wait its full acquisition 
time after the aircraft is again in view. 

5. Terrain Simulation. 

In any computer model of a ground combat environment it is essential 
to have a realistic simulation of the terrain in order to determine the 
mask angles of ground weapon positions. To do this the simulation must 
provide the elevation of any point in the combat area. Present computer 
methods of terrain simulation either use very specialized equipment or 
involve prohibitive amounts of computer storage. The common method is 
to divide the terrain of interest into equal areas and assignan appropriate 
elevation to each. The assigned elevation might be the average elevation 
of the area. In this method it is assumed that all points in a given area, 


for the purpose of simulation, have the same elevation. This is a 


(o>) 





crude approximation of terrain features for large areas. The approximation 
approaches reality as the size of the areas decreases. However, since 
the number of the areas considered is inversely proportional to the size of 
the areas, as the approximation of the terrain approaches reality the com- 
puter storage required becomes very large. 

For the model developed here the terrain features are expressed in 
polynomial form, z=f(x,y). The general problem is that of fitting points 
in three dimensions to a surface represented by the polynomial, z=f(x,y). 

A method of forming this polynomial using the principle of least squares, 
given the set of data points { (x,y,z) } , had been programmed in FORTRAN-60 
for the CDC 1604[2]. This program was converted to FORTRAN-63 and a suc- 
cessful attempt was made in fitting a polynomial to a fairly rough surface 
(mountainous terrain.) The polynomials that best fit mountainous terrain 
whose elevations ranged from 800 to 1300 feet have approximately 400 terms 
with the maximum degree of both x and y near 20. 

The program designed to fit the polynomial to the terrain required the 
input points be from an equally spaced grid in the x,y plane. For this thesis, 
areas of 1000 meters on a side were arbitrarily chosen from maps of actual 
terrain. In order that the fit be acceptable on the edges of the grid square, 
an area of at least 1200 x 1200 meters with the desired 1000 meter Square 
centered inside had to be considered. Consecutive x and y coordinates of 
each individual contour line were recorded and a program was written to 
produce equally spaced elevations from the consecutive data points. How- 
ever, because of the great amount of storage required, these two programs 
could only be incorporated into one by the uSe of overlays. An overlay is 
a feature of FORTRAN-63 designed to increase the use of the core storage 
of the computer by using magnetic tape as auxiliary storage. 

The program to compute z= f(x,y) from the consecutive (x,y) points 


of actual contour lines consists of a main program and three overlays. 





The sole purpose of the main program is to call the individual overlays 
into operation. The first overlay, MAPPING, computes the equally 

spaced elevation points using the (x,y) points from the contour lines as 
input. The second overlay, SURFIT, by using the points computed in the 
first overlay, computes the polynomial z= f(x,y) and records the co- 
efficients on punched cards. The second overlay also computes the ele- 
vations of a grid for use in the third overlay. The third overlay, PLOTTER, 
generates data and stores it on a magnetic tape suitable for a CDC-160 
program used to drive a CalComp 165 graph plotter. The contour map drawn 
can then easily be compared to the original map to estimate the goodness 
of fit of the polynomial. The drawing of this map is optional and is con- 
trolled by an input parameter to the main program. A detailed description 
of each overlay follows. 

Overlay 1 determines the elevation of the terrain at equally spaced 
points of a grid in the x,y plane. The program was written by T. K. Rode- 
burg of the Norwegian Defense Institute for Stanford Research Institute, 
Fort Ord, California. This overlay uses as input the lists of the x and y 
coordinates at successive points on each contour line read from an en- 
larged map of the actual terrain. The following steps are taken by this over- 
lay in the process of determining the elevation of the grid points: 

(1) The point coordinates are corrected for distortion caused by un- 
equal magnification of the original map. 

(2) Rig meals spaced vertical lines are determined and the coordinates 
of the intersection of the contour lines and these vertical lines are calcu- 
lated. 

(3) The elevation at these intersections are calculated by linear 
interpolation. 

(4) By further linear interpolation, elevations at equally spaced 
points on the vertical lines are found, thus obtaining the elevations of 


points on a grid in the x,y plane. 





(5) Using the original data points, the above interpolation procedure 
is repeated at equally spaced horizontal lines. 

(6) The average of the two values thus obtained is used as the ele- 
vation at each point of the grid. 

Overlay 2 fits, in the least squares sense, a polynomial function of 
two independent variables, x and y, to values of a dependent variable, z, 
specified at points on a rectangular grid in the plane of the independent 
variables. That is, given a set of ground elevations, z at grid coordinates 

N M 

fo ges a 


X'y™, which best 


fits the data. The use of orthogonal polynomials leads to a particularly 


x and y, this overlay finds the polynomial, 


simple system of linear equations rather than the ill-conditioned system 
which arises from the usual least Squares normal equations. It also provides 
a measure of the improvements resulting from each new term included which 
further leads, in this algorithm, to an automatic selection of a "best" degree 
polynomial function as determined by Gauss' criterion [1]. It therefore de- 
termines the optimum values of N and M as well as the coefficients. 

This algorithm was programmed in FORTRAN-60 for the CDC 1604 by 
Mrs. Jean Bow of the USNPGS from an ALGOL program [2]. Minor modifi- 
cations were made to convert it to a FORTRAN-63 overlay. The FORTRAN-60 
version is a very general program whereas the overlay has had many of the 
generalities removed for ease of handling inputs which in tum reduced the 
computer running time significantly. A FORTRAN listing of this program 
appears in Appendix IV. 

The following steps are taken by Overlay 2 in the process of determin- 
ing an approximating polynomial: 

(1) the variables are normalized by subtracting off their means, 

(2) the orthogonal polynomials are evaluated , 

(3) the contribution of each orthogonal polynomial to the minimization 


of the residuals is evaluated, 





(4) Gauss' criterion is applied to determine the degree of poly- 
nomial which yields the closest fit to the given data. 

(5) the orthogonal polynomial coefficients are evaluated, 

(6) the dependent variables are evaluated using the approximating 
polynomial, and 

(7) the terrain elevations utilizing a 50 meter grid are evaluated for 
input into Overlay 3. 

Overlay 3 generates lists of x and y coordinates of successive 
points on each contour line and stores these coordinates on an output 
tape in a form suitable for a CDC-160 program uSed to drive a CalComp 165 
plotter. This tape is used as an input tape for the plotter which then recon- 
structs the terrain by connecting the points with straight line segments. 

The FORTRAN-63 program used in Overlay 3 was originally developed 
for use with the IBM 7090 system and an x,y plotter, to give a complete 
traditionally styled contour map. This program was converted to FORTRAN- 
63 by Mr. R. R. Hillery of the USNPGS, for use on the CDC 1604. Subrou- 
tines used to compute the points to be plotted were provided by Dr. M. O. 
Dayhoff of the National Biomedical Research Foundation, Silver Spring, 
Maryland [3]. Subroutines used internally to write the output tape were 
programmed by Mr. J. R. Ward of the USNPGS for use with a CalComp 
165 graph-plotting routine, DRAW. 

Figure 3 is an actual contour map from which approximately 1200 input 
points have been read. A polynomial of 400 terms was fitted and the ele- 
vations of a grid (spacing of 50 meters) were computed to plot Figure 4. 
The coefficients of this polynomial ranged from 904.18 for the constant 
term to -4.593 x 10 . for the aa term. 

An analysis of a table of differences showed the maximum difference 
of 625 actual and computed elevations to be 3] feet in an elevation range 


of 900 to 1200 feet. 
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FIGURE 3 


ACTUAL CONTOUR MAP OF GRID SQUARE 7174 


U. S. ARMY HUNTER-LIGGETT MILITARY 
RESERVAT 
MS 


—e ION 

SS 

oe all 
S 6 \ 


ae 








FIGURE 4 
COMPUTER SIMULATION OF GRID SQUARE 7174 
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Figure 5 shows a representative fit at the junction of two of the grid 
squares approximated by separate polynomials. The accuracy with which 
the grid squares fit together is determined by the accuracy of the input 
Gad. 

6. Analysis of the Model 

PENAIR has been checked out with the use of a simple scenario. A 
3000 x 4000 meter mountainous section of the U. S. Army Military Reser- 
vation was selected for simulation. Three 10-men enemy squads armed 
with M-14 rifles were situated as illustrated in Figure 6. A UH-IB heli- 
copter was assigned the mission of delivering supplies from an airfield 
at position a to an outpost at position b . Two test runs were made with 
altitudes of the helicopter of 1000 feet and 500 feet above the terrain. In 
both tests all 18 evasive maneuvers were checked. Some representative 
flightpaths generated by the program are illustrated in Figures 6, 7, 8, and 
9. The survival probabilities of the aircraft on the various attempts were 
computed and are shown in Tables 1 and 2. All of the computed proba- 
bilities are realistic and are, in a comparative sense, what one would 
expect. 

A sensitivity analysis was then conducted on the rate of climb and des- 
cent and changes in aircraft speeds. A helicopter was simulated flying at 
100 knots with its maximum rate of climb and descent taking the values of 
500,1000 and 2000 feet per minute. The effects of these three rates on 
nap of earth flight are illustrated in Figure 10. It should be noted that 
Figure 10 illustrates how the program generates the flightpath to overshoot 
the crest of a hill before descending just as a pilot tends to do. Figure 10 
also shows, as should be expected, the greater the rate of climb and des- 
cent, the closer the aircraft remains to the earth. Figure 1] illustrates an 
aircraft flown with a maximum rate of climb and descent of 1500 feet per 
minute and speeds of 100, 200 and 300 knots. The resulting flightpaths also 


appear realistic. 


13 


The flightpath computed by PENAIR illustrated in Figure 12 is of a 
helicopter making a level tum at 70, 90 and 110 knots. Here it is noted, 
for the computed flightpath, the radius of turn increases as the airspeed 
increases. 

The last analysis made was the comparison of computed results with 
actual results. Figure 13 is the simulated flightpath and terrain of an 
actual flightpath flown by a UH-1B over an actual terrain. The 44 weapon 
positions illustrated were actually manned by U. S. Army personnel. The 
observed time the helicopter was in view of each of these positions is 
listed in Table 3. The times in view for the simulated terrain and flight- 
path were computed and are also listed in Table 3. It is believed that with 


better terrain data an even closer comparison could be obtained. 
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FIGURE 6 


A STRAIGHT AND LEVEL FLIGHTPATH 
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FIGURE 7 


A FLIGHTPATH OF AN AIRCRAFT CLIMBING STRAIGHT AHEAD 


FIGURE 8 


A PLIGHTPATH OF AN AIRCRAFT TURNING AWAY FROM A FIRING WEAPON 
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FIGURE 3 


A NAP OF THE EARTH FLIGHTPATH 
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Summary of Survival Probability of A UH-1B Helicopter Fired on by Three 
10-Men Squads of M-14 Rifles 


Initial Speed 90 Kts Initial Altitude 650 Meters 


Survival Probabilities 





Evasive Maneuver Evasive 


Commenced after receiving fire Speed WpnA WpnB WpnC Total 





Fly straight and level 70 .8576 .8078 1.0000 .6928 
90 .8723 .8649 1.0000 .7545 

110 .8723 .87238 1.0000) 27el0 

Climb straight ahead 70 .8723 .8723 1.0000 .7610 
90 .8723 .8723 1.0000 .7610 

110 .8723 .8723 1.0000 .7610 

Turn away, level 70 .9026 1.0000 1.0000 -9026 
90 .9026 1.0000 1.0000 -9026 

110 -9026 1.0000 1.0000 .9026 

Turn away and climb 70 .9026 1.0000 1.0000 .9026 
90 .9026 1.0000 1.0000 .9026 

110 .9026 1.0000 1.0000 .9026 

Dive and fly NOE 70 .8359 .7252 1.0000 .6062 
90 .8723 .8031 1.0000 .7005 

110 .8723 .8242 1.0000 .7190 

Turn away and fly NOE 70 :8798 1.0000 1.0000 «8298 
90 .8874 1.0000 1.0000 .8874 

110 .8874 1.0000 1.0000 .8874 

Table ] 


ae) 


Summary of Survival Probabilities of UH-IB Helicopter Fired on by Three 


10-Men Squads Armed with M-l14 Rifles 


Initial Speed 90 Kts 


Evasive Maneuver 
(Commenced after receiving fire) 


Fly straight and level 


Climb straight ahead 


Tur away, level 


Turn away and climb 


Dive and fly NOE 


Turn away and fly NOE 


Evasive 


70 
90 
110 


70 
90 
110 


70 
90 
110 


70 
90 
110 


70 
90 
110 


70 
90 
110 


Table 2 
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Initial Altitude 450 Meters 


Survival 
Speed WpnA WopnB 


~8309 


Bieleilé (i 


.8874 
.8874 
8874 


.8798 
.8798 
.8798 


-/ 556 
O07 
. 3306 


2O796 
-o7 96 
ste feshs: 


Values not determined 
at these speeds 


1766 Ie0Ceu 


.8242 


0000 
0000 
0000 


.0000 
0000 
0000 


v5 01 
.9664 
.Q0000 


-Q0000 
0000 
0000 


HG, 


1 
0000 
0000 


Values not determined 
at these speeds 


0000 


0000 


-0000 
-0000 
.0000 


.0000 
0000 
.0000 


0000 
.0000 
.0000 


Probabilities 
WpnC Total 
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.8874 
8874 
.8874 
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FIGURE 10 


COMPARISON OF RATES OF CLIMBS AND DESCENTS IN A NAP 
OF THE EARTH FLIGHT AT A CONSTANT SPEED OF 100 KTS. 
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COMPARISON OF AIRCRAFT SPEEDS ON A NAP OF THE EARTH FLIGHT 
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FIGURE 12 


COMPARISON OF THE AIRCRAFT'S RADIUS OF TURN 
AT 3 AIRSPEEDS, 70, 90 and 110 KTS. 
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FIGURE 13 


COMPUTER SIMULATION OF AN ACTUAL FLIGHTPATH 
FLOWN BY A UH-1B HELICOPTER 
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Comparison of Actual and Computed Time Intervals in View of 44 Weapon Sites 


Time Interval = .6 Sec. 


Wpn. Time in View Wpn Time In View Wpn Time In View 
No Actual Computed No Actual Computed No Actual Computed 


1 0 0 31 0 0 
2 0 9 32 22 39 
3 30 30 33 0 0 
4 0 34 89 88 
5 0 35 93 91 
6 0 59 36 69 76 
7 0 0 a7 66 58 
8 0 40 38 49 64 
g 0 0 39 16 36 

10 0 17, 40 0 16 

al 9 2 4} 8 8 

12 0 0 42 56 56 

13 26 39 43 19 15 

14 40 63 44 103 102 

15 55 73 





Table 3 
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7. Conclusions and recommendations. 

The simulation of the terrain obtained by this program has proven to be 
extremely realistic. (See Figure 4, Section 5, page 1] ). The junctions of 
Squares approximated by separate polynomials are as good as the data used 
for the approximation. (See Figure 5, Section 5, page 12). This program was 
written in its present form because of the availability of input data to the 
author in the format described in Appendix V. If the actual elevations could 
be obtained using an equally spaced grid for the x and y coordinates, the 
first overlay could be eliminated. It is felt that this is feasible and would 
make the preparation of the input cards much simpler. 

Among the uses of the simulation of the terrain by this method include 
ground to ground combat and ground to air combat environments. Maskangles, 
terrain profiles and actual distances between points in a mountainous terrain 
are all easily determined. The uses of the PENAIR program include the deter- 
mination of optimal flightpaths across hostile terrain, the survival proba- 
bility and optimal maneuvers to perform if the aircraft is fired upon. 

The results in Section 6 all tend to illustrate the realism of the aircraft 
penetration model. Whether the actual kill probabilities used are exact is 
not extremely important since an appropriate defensive maneuver for an un- 
armed aircraft to initiate upon receiving ground fire may be obtained by ob- 
serving the relative survival probabilities of several maneuvers. 

It is believed that with additional programming this model could be ex- 
tended to include formations of two or more aircraft, thus greatly increasing 
its usefulness. In addition, since the terrain is simulated by mathematical 
functions and since the line of sight between aircraft and weapon may be 
described by a linear relationship, it appears possible to solve the two rela- 
tions for a point of intersection. Upon accomplishing this programming, it 


is believed the computer running time would be reduced significantly. 
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APPENDIX |] 


Flow Charts of the PENAIR Program 


In the flow charts of this appendix, the following conventions 


are used: 


saad 


41.8 420 


©) a 


This symbol is used for a READ 


statement. 


This symbol is used for a PRINT 


statement. 


A rectangle indicates one or more 
FORTRAN statements. A number such 
as 4.1.8 indicates the statement 

in this rectangle is further explained 
by Flow Chart 4.1.8 in this appendix. 
A number such as 420 is the statement 
number as it appears in the FORTRAN 
listing in Appendix II. 


This symbol is used for an IF statement. 


This symbol indicates the logic of this 
Flow Chart is continued at the beginning 


of Flow Chart 5.1 of this appendix. 
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1 OO 


CONT/NUE. 


1 


L+/,/2 


aN 
S 
= 


This symbol indicates a connector between 
two parts of the same Flow Chart of this 


appendix. 


This symbol indicates the logic of this 
Flow Chart is continued where the same 
symbol appears in Flow Chart 4.1 of this 
appendix. 


This is the symbol for a CONTINUE state- 


ment whose statement number is 10. 


This is the symbol for a computed 
GO TO statement. 


This symbol is used for the statement 


DO 10 I1=1, l2 


This symbol indicates the beginning or 


entry point for the Flow Chart. 


This symbol indicates the end or exit 


point for the Flow Chart. 
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APPENDIX II 


FORTRAN Listing of the PENAIR Program 


PENAIR consists of the main program and four subroutines. The 
variable names used in the PENAIR program are defined in the comments 
appearing in the FORTRAN listing. The variable names of input and out- 


put parameters are explained in Appendix III. 
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Appendix III 
Inputs and Outputs of PENAIR 


Data cards for the PENAIR program should be prepared in the following manner 


Data Number Variable Columns Description 
Group ene ards Name of Card of Variable 
I l NRUN 1-10 A right justified integer, 


the number of runs to be 


performed. 


2. l KMAX L=10 A right justified integer, 
the number of squares used 
to approximate the terrain, 
KMAX < 12. 

XMIN 11-20 A right justified floating 
point number, the X co- 
ordinate of the left edge of 
the terrain. All coordinates 
are in meters, that is the 
distance between the two 
xX coordinates 72455 and 
72454 is one meter. 

YMIN 21-30 A right justified floating 
point number, the y co- 
ordinate of the bottom edge 
of the terrain. 

XMAX ol—40 The x coordinate of the right 


edge of the terrain. 
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Variable Columns 
Name of Card 
YMAX 41-50 
thet 51-60 
TPD 61-70 
Xx 1=10 
11-20 

Z 21-30 
X2 31-40 
YZ 41-50 
GMAX 1-10 
RC 11-20 
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Description 
of Variable 


The y coordinate of the 
top edge of the terrain. 

A right justified floating 
point number, the time 
interval of the game 
(seconds), TI>0O. 

A right justified floating 
point number, the dimen- 
sion of the map square 
approximated by one poly- 


nomial (meters) . 


Right justified floating 

point numbers, the begin- 
ning coordinates of the 
aircraft. If IMARK = 3,4, or 
6, y must equal YMIN. 

The x and y coordinates of 
the desired destination of 
the aircraft. If IMARK = 7 


or 8, this card may beblank. 


A right justified floating 
point number, the maximum 
positive G force to be used 
by the aircraft. 

A right justified floating 
point number, the aircraft's 
maximum rate of climb 


(feet/minute) 





Data Number Variable Columns 


Group of Cards Name of Card 
RD 71=30 

GS 31-40 

5 i IMAX ae 

6 l WPNAWPNB 1-16 
JMAX Z=25 

TA 26-30 

ET Soleo 
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Description 
of Variable 


A right justified floating 
point number, the aircraft's 
maximum rate of descent 
(feet/minute) . 

A right justified floating 
point number, the original 


airspeed of the aircraft. 


A right justified integer, 
the number of types of 


weapons, IMAX< 10. . 


The name of the weapon, 
e.g., 90 Cal. MG 

A right justified integer, 
the number of weapons of 
this type, JMAX < 20. 

A right justified integer, 
the time required to acquire 
the target by this weapon 
(seconds), TA20. 

A right justified integer, 
the firing time of this 
weapon before reloading 
is required (seconds) , 


File 6). 





Data Number Variable Columns Description 
Group of Cards Name of Card of Variable 
RLT 36-40 A right justified integer, 
the time to reload this 
weapon (seconds) RLT 20. 
RMAX 41-50 A right justified floating 
point number, the maximum 
range of this weapon(meters) 
RMAX>0O. 
RMIN 51-60 A right justified floating 
point number, the minimum 


range of this weapon(meters) . 


RMIN 20. 
7 2 P 1-8 The probability of kill by 
9-16 this weapon in one time 


interval at ranges of 


7280 1g@, 200, 300.,.....9eW00 


meters. 
1-8 The second card of this 
9-16 group contains the kill 


probabilities at ranges of 
72-80 1100 to 2000 meters in 100 


meter intervals. 


8 JMAX NW 1-2 The number assigned to the 
weapon, 1<sNW<s99. 


XW 11-20 The x, y and z coordinates 
YW 21-30 of this weapon, ZW may be 
ZW 31-40 blank since PENAIR will com- 


pute the weapon's elevation. 
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NOTE: Data groups 6, 7 and 8 should be repeated IMAX times. 


Data Number Variable 
Group of Cards Name 
9 il MMAX 
NMAX 
10 As many cards PHI 


as necessary 


Columns 
of Card 


Lstd 


11=20 


Description 
of Variable 
A right justified integer, 
the degree of x in the poly- 
nomial that approximates 
the first map square. 
A right justified integer, 
the degree of y in the 
polynomial that approximates 


the first map Square 


The coefficients of the 
first polynomial appear on 
this and the following cards 
of this data group ina 
4E20.10 field as obtained 
from the terrain simulation 


program. 


NOTE: Data groups 9 and 10 should be repeated KMAX times. 


11 4 AFIELD 


lee: 


These four cards contain 

the four lines of the title 

to be printed at the top of 
the output. 





Data 
Group 


12 


Number Variable Columns Description 
of Cards Name of Cards of Variable 
l IMARK 10 The-value of this flag, deter- 
mines the flightpath to be flown. 
IMARK = 1, upon receiving fire, 
fly straight and level, 
= 2, climb straight ahead, 
= 3, turn away, level, 


= 4, turnawayandclimb, 


| 


= 5, dive and flynap of 
the earth, 

= 6, turn away and fly 
nap of the earth, 

= 7, the flight path is 
completely pre- 
determined andread 
in as data groupl3, 

= 8, the x and y coordi- 
nates are read in 
as data group 13, 
and a napofthe earth 
flightpath is computed. 

IPUNCH 20 IPUNCH = 1, the flightpath will 

will be punched on 


cards. 
= 0, the flightpathwill 


not be punched. 
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Data Number Variable Column 
Group of Cards Name of Card 


GSNEW 21-30 


Description 
of Variable 


A right justified floating 
point number, the ground 
speed to be attained by the 


aircraft after receiving fire. 


NOTE: Data groups 1] and 12 will be repeated NRUN times. 


13 As many NPATH 6-7 
cards as 
necessary 
but not 
greater XT 18-26 
than 300. VT 27-35 
ZH 36-44 


Any positive integer, used 


to identify the flightpath, 


Right justified floating 
point numbers, each of the 
cards of this group contain 
xX, y, and z coordinates of 
the aircraft in consecutive 


time intervals. 


NOTE: A blank card should be inserted at the end of this data group. 
This data group will be used only if IMARK = 7 or 8. If 
IMARK = 8, ZH may be blank. Data group 13 should follow 


its corresponding data groups 1] and 12. 
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Output of the PENAIR Program 


A sample output from the PENAIR program is illustrated in Figures 
14 and 15. The input parameters of the weapons are printed with the 
weapon positions and elevations. The second page of the output (Figure 
15) is a summary of time in view and aircraft survival probability. The 
time in view listed here is in number of time intervals in view. If the 
time interval used in this example had been 0.5 second, the aircraft 
was in view of weapon number 49, located at (70080, 76940, 339), for 
28 seconds (56 x 0.5). The survival probability of the aircraft from this 


weapon was .999]l. 
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PIRST PAGE OF OUTPUT FROM THE PENAIR PROGRAM 
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APPENDIX IV 
FORTRAN Listing of the TERRAIN Program 


TERRAIN consists of a main program and three overlays. The control 
cards for the COOP-Monitor of the CDC 1604 are included in this listing 
and are explained in Appendix V. The variable names used in the TERRAIN 
program are defined in the comments appearing in the FORTRAN listing. 
The variable names of input and output parameters are explained in Appendix 


V. 
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APPENDIX V 
Inputs and Outputs of the Terrain Simulation Program 


The following information is presented for potential users of this pro- 
gram: 


A. Tape units 


(1) The COOP control cards are shown in the FORTRAN listing (See 
Appendix IV.) A dash (-) in column one indicates a 7,9 punch, while an 


open parenthesis (() indicates an 11,0,7,9 punch. 


(2) Logical tape unit 49 is the graph tape and 45 refers to the type- 


writer. 


(3) Scratch tapes 56 and 57 are necessary and therefore a total of 


eight tape units are required. 
B. Input Formats 


The preparation of the input data cards is explained on the following pages 


of this appendix. 


Il ath 





Data Number Variable Columns 
-Group of Cards Name of Card 
| 2 LTEILE 1-48 
Z | NL 1-95 
3 NL IT j-5 
LABL(ID) 7-10 
4 ] NMAX 1-9 
MMAX 6-10 


NEWCTL T¥alo 


le 


Description 

of Variable 
These two cards contain the 
two line title desired for 
the first page of the output 
and the bottom of the plotted 


map. 


A right justified integer, 
the number of contour lines 


requiring a label. 


A right justified integer, 
the number of the contour 
level to be labeled, 
l<eIIs 30 

A left justified integer, 
the elevation of contour 


level II. 


A right justified integer, 

the maximum degree of x to 

be considered 

A right justified integer, 

the maximum degree of y to 

be considered. 

A right justified integer, 

NEWCTL =0, no contour map 
is desired. 

=l1, an outputtape 

will be generated 
for the CalComp 
165 plotter. 





Data Number Variable 


Group of Cards Name 
se) 1 S 
A 
6 il XTMIN 
XTMAX 
CND 
7 iL GX 


Columns 
of Cards 


Lalo 


LJ=20 


Ll=20 


21-30 


le2 
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Description 

of Variable 
A right justified floating 
point number, the spacing 
desired between grid lines. 
S must be chosen so that 
neither the number of vertical 
nor horizontal grid lines 
exceed 45. 
A right justified floating 
point number, the linear 


dimension of the map Square. 


A right justified floating 
point number, the elevation 
of the lowest contour in 

the map square. 

A right justified floating 
point number, the maximum 
elevation in the map Square. 
A right justified floating 
point number, the contour 
interval. The first contour 


is the value C = XTMIN. 


The first two digits of a 
five digit floating point 
number, GX * L000 = the x 
coordinate of the lower 


left corner of the map square. 





Data 
Group 


Number 
of Cards 


Variable 
Name 


Gr 


XP (1) 


Yeu) 


XP(2) 


YF(2) 


XP (3) 


XP( 3) 


XP(4) 


XP(4) 


Columns 
of Cards 


3-4 


Neh} 


16-20 


ei-2o 


26-30 


31-35 


36-40 


41-45 


46-50 
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Description 

of Variation 
The first two digits of a five 
digit floating point number, 
GY * 1000 = the y coordinate 
of the lower left comer of the 
map square. 
A right justified floating point 
number, the x coordinate of 
the lower left comer of the map 
square measured on an arbitrary 
scale. XP (1) usually should 
equal 0000. 
The y coordinate of the lower left 
comer measured on an arbitrary 
scale. YP(l) usually should equal 
0000. 
The x coordinate of the lower 
right corner, 
The y coordinate of the lower 
right corner. 
The x coordinate of the upper 
right comer. 
The y coordinate of the upper 
right corner. 
The x coordinate of the upper 
left comer. 
The y coordinate of the upper 


left corner. 





Data Number 
Group of Cards 
8 As many 
cards as 
necessary 


Variable | Columns 
Name of Card 

1-4 

FZ 5-8 

NC 9-)0 

X(1) 1l-15 

Y¥() 16-20 
21-80 


Description 
of Variable 


Any identifying ccde desired 
Col. 1-4 are not read. 

A right justified floating 
point number, FZ +*« 10 = 

the elevation (feet) whose 
coordinates appear on this 
card. 

A sequence number starting 
at 01 and increasing by one 
on each card following whose 
coordinates are thoseofthe 
Same contour line. 

A right justified floating 
point number, the x coordi- 
nate of a point on the contour 
line. 

A right justified floating 
point number, the coordinate 


of a point on the contour line. 


F 5.0 fields of additional 
points on the contour line 
in the same manner as 


Col. 11-20. 


NOTE: One blank card is inserted after this group of data cards (eighth 


data group) for each contour line. One additional blank card is 


inserted if a new seventh data group is read. 
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Data 
Group 


NOTE: 


Number Variable 
of Cards Name 


Columns 
of Card 


Pm 


Description 

of Variable 
The figure 999 should appear 
in these colums to Signify 
that this is the last card 
of the data deck. 


There should be no blank cards between the last efghth data group 


and the ninth data group. 
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Example: 

It is desired that grid square 7175 of the U. S. Army Hunter-Liggett 
Military Reservation be approximated by a polynomial. Figure 17 illus- 
trates the data cards necessary 

The words 

HLMR CLOSE TERRAIN 

GRID SQUARE 7175 
will appear at the top of the first page of the printed output and at the 
bottom of the plotted map. The minimum and maximum elevation appear- 
ing in this square are 900 and 1200 feet respectively. The contour interval 
is 20 feet. The maximum degree of X and Y to be considered is 20. The 
grid interval will be 50 meters and the square is 1000 meters on a side. 

Figure 16 illustrates how the original map was enlarged to read the 


contour coordinates. 


(001 2, 3485) (3455, 3440) 


(0000, 0000) (3452, 0008) 
FIGURE 16 


DISTORTION CAUSED BY PROJECTION 
OF A MAP GRID SQUARE 


Distortion caused by projection is apparent in Figure 16 and is cor- 
rected by the first overlay. The corner coordinates are arbitrary and all 
contour coordinates are relative to them. If the points are measured by 


hand, the actual corner coordinates may be used and in this case would 
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FIGURE 17 


A SAMPLE DATA DECK FOR THE TERRAIN PROGRAM 
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be (71000, 75000), (72000, 75000) ,(72000, 76000), (71000, 76000). The 
contour coordinates would then be the actual terrain coordinates. 
In this example the grid square was reprojected a second time with new 


corner coordinates of (0000, 0000), (3444,-0020), (3441, 3467), (-0004,3475). 


Terrain has three types of output, a) printed, b) punched cards, and 
c) a plotted map (See Figure 18). The printed output consists of 1) the x and 
y coordinates of the grid to be used, 2) the input elevations at points on this 
grid as interpolated by Overlay 1,3) the elevations of the same points as 
computed using the approximating polynomial, 4) a table of differences be- 
tween the actual and computed elevations of the grid, 5) the computed ele- 
vations to be used in Overlay 3, 6) the coefficients of the approximating 
polynomial, and 7) the actual mean elevation squared, computed mean 
elevation squared and the maximum differences between the actual and com- 
puted elevations. Although accuracy is not determined in general, if the 
actual mean elevation squared and computed mean elevation squared are 
equal, then computation is exact. In practice they will not be equal due 
to the imprecise nature of calculation. A wide discrepancy indicates ex- 
cessive errors in calculation. 

The first card of the punched card output contains the degree of x and 
y ina 2110 format. The remaining cards contain the coefficients of the 


approximating polynomial in a 4E20.10 format. 
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HLMR CLOSE TERRAIN 
CRID SQUARE 7174 


FIGURE 18 


A SAMPLE OF A PLOTTED MAP 
FOR THE TERRAIN PROGRAM 
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